Search for W and Z Bosons in the Reaction pp 2 jets + 7 at ^/s = 1.8 TeV 
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We present a study of the dijet invariant mass distribution for the reaction pp — > 2 jets+7 + X, 
at a center of mass energy of 1.8 TeV, using data collected by the CDF experiment. We compare 
the data to predictions for the production of a photon with two jets, together with the resonant 
processes pp W/Z + 'y + X , in which the W and Z bosons decay hadronically. A fit is made to the 
dijet invariant mass distribution combining the non-resonant background and resonant processes. 
We use the result to establish a limit for the inclusive production cross section of W/Z + 7 with 
hadronic decay of the W and Z bosons. 



PACS numbers: 

INTRODUCTION 

The production of final states containing W oi Z 
bosons in association with a photon (7) in pp collisions 
at y/s = 1.8 TeV has been studied by the CDF and D0 
collaborations using event samples in which the W and 
Z bosons decay to leptons . Identification of Wj and 
Z"f events in which the W/Z decay hadronically is ex- 
perimentally difficult because of large background from 
direct production of 2 jet -f 7 and three-jet events. A 
study of the general kinematic characteristics of 2 jet -|- 
7 production has been conducted using 16 pb^^ of data 
from CDF Run la, however the dijet invariant mass dis- 
tribution was not investigated 0. 

We present a search for evidence of 7 -f W/Z — > qq fi- 
nal states in 90 pb~^ oipp data. The data were collected 
with the CDF detector during Run lb of the Tevatron. 
A search for 2 jet -I- 7 events was conducted within the 
data subset in which a photon candidate had electromag- 
netic transverse energy greater than 23 GeV. The dijet 
invariant mass distribution was fit to a mixture of boson 
resonance decay and QCD background. 

The methods established in this analysis Q could 
prove useful for identifying similar signals coming from 
X jj decay which are embedded in large QCD back- 
ground. 



THE DETECTOR 

A comprehensive description of the Collider Detec- 
tor at Fermilab (CDF) may be found in 4]. We used 
a coordinate system with z along the proton beam, 
azimuthal angle 0, polar angle 9, and pseudo-rapidity 
rj = —\n tan{9/2). The transverse energy of a particle 
(e.g. photon, electron, jet) was defined as Et = Esm6. 
The primary components relevant to this analysis were 
those that measure jet energies and positions, photon en- 
ergies, positions and profiles, and those that establish the 
pp interaction vertex. 

The central tracking chamber (CTC) and vertex track- 
ing chamber (VTX) were used to measure momenta and 
directions of the charged tracks associated with jets. The 
tracking chambers were located within 1.4 T axial mag- 
netic field. The CTC was a drift chamber which provides 
space point information used to construct the trajecto- 



ries of charged particles. It covered a rapidity range of 
\ri\ < 1.1. The VTX was a time projection chamber po- 
sitioned between the beam pipe and the CTC that pro- 
vided improved interaction vertex measurement with the 
extrapolation of tracks reconstructed in the CTC. 

Scintillator-based electromagnetic (CEM) and 
hadronic (CHA) calorimeters in the central region 
(1^1 5: 1-1) were arranged in projective towers of size 
Af] X w 0.1 X 0.26. The end- wall hadronic calorimeter 
(WHA) and the end-plug electromagnetic (PEM) and 
hadronic (PHA) calorimeters extended the rapidity 
coverage out to \r]\ < 2.4. 

Two additional detector elements were used for photon 
identification. The central strip chambers (CES) were 
multi-wire proportional chambers with segmented cath- 
ode strips. The CES was positioned near the shower 
maximum of the central electromagnetic detector. The 
anode wires measured (p and the cathode strips mea- 
sure 77 for showers in the CEM. The central preradiator 
(CPR) was a set of multi-wire proportional chambers, 
positioned between the magnet solenoid and the CEM. 
It was used to measure the electromagnetic shower pulse 
heights of electron-positron pairs from photons convert- 
ing in a solenoid of thickness 1.1 radiation lengths. The 
CES and the CPR detector systems provided discrimi- 
nation between single photons and multi-photon showers 
arising from tt" and 77 decays. 

EVENT SELECTION 

Events with photon candidates were selected using a 
trigger which required a high threshold of energy de- 
posited in the central electromagnetic calorimeter. Pho- 
ton clusters were reconstructed by combining the energy 
from neighboring CEM cells with a seed cell having en- 
ergy above a threshold of 3 GeV. A photon candidate 
was required to have clustered energy Et > 23 GeV. 

The trigger acceptance for photon candidates has been 
measured as a function of the photon transverse momen- 
tum 1^. The acceptance plateaus for the photons with 
Et > 30 GeV at 0.970 ± 0.006. 

In the offline data analysis, the measured photon en- 
ergies were corrected using an algorithm taking into ac- 
count variations of CEM response within a cell, cell-to- 
cell energy sharing, and a global transverse energy scale 
. A photon candidate was required to be isolated, with 
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less than 15% additional energy {Eqem + Echa) within 
a cone of radius Ai? = + A772 < 0.4 centered 

on the photon direction. In order to reject electrons, a 
tracking isolation requirement was imposed by summing 
the transverse momentum of tracks within the Ai? < 0.4 
cone around the photon direction. The sum of pT of the 
tracks was required to be less than 2 GeV/c. 

The ratio of the energy deposited in the hadronic 
calorimeter {Eh ad) to the energy deposited in the elec- 
tromagnetic calorimeter (Ecem) for a photon candidate 
passing the above cuts must satisfy the requirement 



Ehad 
Ecem 



< 0.055 + 0.00045£;(7), 



with E{'j) in GeV. Typically, for a photon, the above ra- 
tio is less than 10%. A lateral sharing parameter for GEM 
clusters measures the spread of energy over calorimeter 
cells adjacent to the seed cell and selects photons based 
upon the spread expected from measurements in an elec- 
tron test beam. 

The photons had to pass fiducial cuts that require suf- 
ficient shower containment in the GES chambers. The 
highest energy strip cluster and the highest energy wire 
cluster were chosen to estimate the position and trans- 
verse profile of a photon candidate. Any additional clus- 
ters in the GES within the boundaries of the calorimeter 
energy cluster had to have energy less than 2.39 GeV 
+0.01 Et{i) GeV. 

Multi-photon backgrounds from tt'^ and r] decays are 
suppressed by requiring the transverse profile of the 
shower energy observed in the GES to match to a ref- 
erence profile measured using an electron test beam. A 
partial separation of direct photons from multi-photon 
backgrounds was made based on the quality of the shower 
shape agreement. Additional details of the photon cuts 
may be found in 0. 

The energy from identified photons is removed from 
the calorimeter cells and the remaining energy is clus- 
tered to reconstruct jets. Hadronic jets are identified us- 
ing a cone clustering algorithm The energies of jets 
clustered with a cone size of Ai? < 0.7 were corrected for 
calorimeter nonlinearity and cell-to-cell non-uniformity. 
A correction takes into account energy deposition from 
particles not associated with the parent parton falling 
within the jet clustering cone, and an out-of-cone correc- 
tion for the underlying event accounts for energy from 
the jet which is radiated outside the jet clustering cone. 

The sample of events with exactly two jets having Et 
greater than 15 GeV and \rj\ < 1.1 is used. The restricted 
rapidity range insures that the reconstruction efficiency 
for tracks associated with these jets is uniform. Any ad- 
ditional jets present in the event within \ri\ < 2.4 must 
have Et < 10 GeV. Each of the two central jets must in- 
clude at least one well reconstructed track, which is used 
to calculate the vertex position. The vertices determined 
from the tracking information from the two jets in an 



event must lie within 10 cm of each other. This cut re- 
moves instances of jets from different pp cohisions during 
the same beam-crossing. The interaction vertex is taken 
to be the average of the two jet vertices and must he 
within 60 cm of the center of the detector. This vertex, 
together with the information from the GES, is used to 
define the photon trajectory. 

A total of approximately 2.7 million triggered events 
was recorded. A sample of 9493 events satisfied all pho- 
ton and jet selection criteria. The photon Et distri- 
bution is shown in Figure 1. Figure 2 shows the Et 
and 77 distributions for the two leading jets in the event, 
ordered by Et- The first jet has average Et of 39 
GeV, while the second jet has average Et of 23 GeV. 
A plot of the dijet invariant mass distribution is shown 
in Figure 3. The dijet invariant mass mjj is defined as 



{El + E2Y - {Pi -t- where E., and Pi are the en- 
ergies and momenta of the two leading jets. 




50 75 100 

Photon E-r (GeV) 



FIG. 1: Photon Et Distribution from 7 -I- 2 jet events. 



EVENT GENERATION AND SIMULATION 

There are several physical processes which produce the 
2 jet -|- photon final state. Our primary goal is the study 
of resonant W/Z + 7 production with W/Z decaying to 
jets. In addition to these signal events, we also generate 
two QGD backgrounds which contribute to the two-jet 
+ photon events. The first contributor is non-resonant 
two-jet -|- photon production and the second is three- 
jet production where one of the jets is misidentified as a 
photon. 
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to the data. 
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FIG. 2: The Et and rj distributions for the two leading jets. 




Dijet Invariant Mass (GeV/c ) 



FIG. 3: Dijet mass spectrum from 2 jet +7 data. There 
are 2656 events in the rrijj region from 60 to 126 GeV/c^. 
The solid fit corresponds to the background + signal fit as 
described in Section 6. 



A CDF detector simulation produces event records 
with the same structure as the data. These generated 
events are then subject to the same photon and jet selec- 
tion cuts, and detector geometric acceptance, as applied 



Signal Simulation 

The signal processes — > W + j and pp Z + 
in which the W/Z decay to 2 jets, are modeled us- 
ing a leading order (LO) matrix element calculation 
[Sl which includes contributions from initial and final 
state inner bremsstrahlung processes. The Monte Carlo 
simulation produces photons with Et > 15 GeV and 
AR between the photon and partons greater than 0.4. 
We also require that the invariant mass of the quark- 
antiquark pair from the W/Z decay is greater than the 
boson mass minus three times the boson natural width 
{VS > My — SFy). This latter cut tends to suppress final 
state bremsstrahlung and produces a characteristic W/Z 
resonant mass peak that simulates the signal searched for 
in our data (see Section 6). The factorization scale was 
set equal to the square of the colliding partons center of 
mass energy The MRSA' parton distribution func- 
tions were used. The parton-level events are run through 
HERWIG for parton shower evolution and hadronization 
[inj |. The resulting dijet mass distributions from W and 
Z decay, together with the combined distribution, are 
shown in Figure 4. 
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FIG. 4: The invariant mass distributions of the two leading 
jets from the LO and Zj simulations, together with their 
sum (solid line). The combined invariant mass was fit with a 
double Gaussian as described in Section 5. 



The W7 and Zj signal acceptance is calculated with 
a full CDF detector simulation. The details of the ac- 
ceptance calculations are presented in Tables 1 and 2 for 
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the and Z7 events respectively. The uncertainty in 
the acceptances for each cut were estimated by varying 
the upper and lower limits on the cut parameters and 
measuring the effects on the number of events passing 
the cuts. The errors in the Table 1 and Table 2 represent 
cumulative errors on the acceptance, where all the errors 
of the subsequent cuts are added in quadrature. 



TABLE I: W7 event selection acceptances 



W7 sample 


events 


% 


error % 


generated 


97298 


100 




photon I77I 


55487 


57.03 


+ 2.40, -2.83 


photon ID 


25939 


26.67 


+ 1.13, -2.02 


photon trigger 


9567 


9.83 


+ 0.62, -0.83 


jet Et + \ri\ 


1691 


1.74 


+ 0.30, -0.29 


extra jet 


1691 


1.74 


+ 0.31, -0.30 


tracks 


1677 


1.72 


+ 0.31, -0.30 



TABLE II: Z7 event selection acceptances 



Z7 sample 


events 


% 


error % 


generated 


102210 


100 




photon \ri\ 


68563 


67.08 


+ 2.55, -3.16 


photon ID 


37080 


36.28 


+ 1.38, -2.78 


photon trigger 


15294 


14.96 


+ 0.86, -1.45 


jet Et + 77 


2959 


2.90 


+ 0.47, -0.49 


extra jet 


2959 


2.90 


+ 0.48, -0.51 


tracks 


2946 


2.88 


+ 0.48, -0.51 



The photon 77 cut was varied from 1.0 to 1.2. The 
photon ID cut combines the number of charged tracks 
associated with a photon cluster, the sum of transverse 
momenta px of the charged tracks, and the photon isola- 
tion. The number of charged tracks allowed in a photon 
cluster cone was varied from zero to two. The effect of 
this change was propagated through to the photon isola- 
tion cuts, and the resulting variation in the acceptance 
was used as an estimate of the systematic uncertainty 
in the photon identification. Finally, we estimated the 
uncertainty in the acceptance due to the photon trigger. 
The key threshold points in the turn-on function ^ were 
varied by ±5% from the nominal value. For example, 
the photons with 23 GeV < Et < 26 GeV originally had 
22% probability to pass the trigger. The high acceptance 
cut lowers this probability to 17%, while the low accep- 
tance cut raises this probability to 27%. In addition, the 
photons with 26 GeV < Et < 30 GeV had 77% proba- 
bility to pass the trigger. Finally, photons with Et > 30 
GeV were passing the trigger with 97% probability. The 
resulting effect on the acceptance was taken as a system- 
atic uncertainty. Tables 1 and 2 summarize the photon 
selection acceptances and their uncertainties. 

The systematic uncertainties due to jet selection crite- 
ria can be divided into three categories: those based upon 



the cuts used to select the two high Et jets {Et > 15 
GeV, and < 1.1); those caused by rejection of events 
having additional low Et jets ( additional jet can not 
have Et > 10 GeV if within \t]\ < 2.4); and those due to 
the use of charged tracks in the jets to define the event 
vertex. The uncertainty caused by the high Et jet se- 
lection was estimated by varying the Et cut from 14 to 
16 GeV and the \ri\ cut from 1.0 to 1.2 (see "jet Et and 
T]" rows in Tables 1 and 2). The effect of the rejection of 
additional jets was evaluated by varying Et from 9.5 to 
10.5 GeV and \ri\ from 2.3 to 2.4 ("extra jet" rows in Ta- 
bles 1 and 2). The final uncertainty on the jet selection 
comes from the use of charged tracks to determine the 
event vertex. The upper number on the tracking uncer- 
tanty was estimated by demanding that each of the jets 
has at least 2 charged tracks which in addition have very 
small distance from the calculated jet vertex. Variations 
in the track selection criteria cause almost negligable un- 
certainty of event acceptance shown in the last rows of 
Tables 1 and 2. The final overall event selection accep- 
tance for W-f events is 0.017 ± 0.003 and for Z7 events 
0.029 ± 0.005. 



Background Simulation 

We use a tree level calculation for the 7 + 2 jets 
background process which includes both prompt photon 
and bremsstrahlung contributions , followed by HER- 
WIG for parton evolution and fragmentation. In addi- 
tion, three-jet events are generated using PYTHIA [l^ . 
with JETSET performing the parton evolution and frag- 
mentation. Each of the three leading jets were tested for 
misidentification as photons. This procedure was based 
on the probability distribution measured from jet data 
which describes how often a jet with a particular Et is 
misidentified as a photon. The data show that a jet with 
Et > 23 GeV has a maximum probability of 8 x 10~^ to 
fake a photon. 

The di-jet mass distribution from the data favors a 
background mixture of 60% two-jet -l-photon and 40% 
three-jet in which one jet fakes a photon. Figure 6 
shows the dijet mass distribution for this combination 
of backgrounds. For masses between 60 GeV/c^ and 126 
GeV/c^ the background distribution is fit well to the form 
dN/drrijj = Aexp^—bmjj) as shown by the solid line in 
Figure 5. This slope is known to be effected by the inclu- 
sion of higher-order QCD contributions, so we leave the 
exponential slope of the background as a free parameter 
in the fit to the experimental dijet mass mass spectrum 
to a combination of signal plus background (see Section 
6 for details). 
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FIG. 5: Combined background of two-jet + 7 events and 
three-jet events in which one jet is identified as a photon. 
The distributions are mixed in the ratio 60% to 40%. The 
exponential fit is presented with a solid line. 



STANDARD MODEL PREDICTIONS FOR VK7, 
CROSS SECTIONS 



Four parton-level event samples were generated using 
the procedure described in Section 4.1: Ty+7 (where 
— > ud), (where W~ du), Zj (where 

Z uu), and Z"f (where Z — > dd). Partons are treated 
as massless in the model. In order to include the W/Z 
decays to heavy quarks, the initial "light" cross sec- 
tions are multiplied by appropriate factors. For exam- 
ple, ud is multiphed by 2.1 to account for cs, cd, 
us and cb decays. The same is valid for the W~ du. 
The cross section for Z ^ uu is multiplied by 2.0 to ac- 
count for Z — > cc, while the cross section for Z — > is 
multiplied by factor 3.0 to include Z ss, bb decays. 
In addition, the final cross sections are multiplied by a 
factor of 1.3 0| to take into account higher order QCD 
corrections. 

The resulting cross section for pp — > + 7 with 
decay to quark-antiquark is 6.35 pb, and for pp 
Z + J with the Z decay to quark-antiquark is 6.52 pb. 
These values pertain to events generated with £'t(7) > 
15 GeV, AR between the photon and each of the quarks 
Ai?(7 — q) > 0.4 and the quark-antiquark mass greater 
than the boson mass minus three times the boson natural 
width. 

With an integrated luminosity of 90 pb~^ we would 
expect Nw-y = 90 x 6.35 x 0.017 = 10^2 events, and 



Nzj = 90 X 6.52 X 0.029 = 171^ events. The total is 27I5 
events in the two-jet + photon data sample in the dijet 
invariant mass region between 60 and 126 GeV/c^. The 
errors on the total number of events were added directly, 
as they result from the same detector systematics in the 
Monte Carlo simulation. The generated W^j and Zj 
events are combined in proportion to their cross sections. 

The dijet mass distribution from the combined W/Z^ 
generated events shown Figure 4 is normalized to a total 
area of one and fit to a double Gaussian 



1 dN 
N dm a 



S{mjj) 



S{mjj) = Sie (="1)" + 826 (2-2)^ 
The fit parameters are = 0.028 (GeV/c^)- 



nil 



and S2 = 0.002 
41 GeV/c2. These 



89.6 GeV/c^ cti = 12.5 GcV/c^ 
(GeV/c2)-i, m2 = 76.3 GeV/c^, aa 
parameters are used in fitting the dijet mass distribution 
from the data to a combination of an exponential back- 
ground and the signal distributions. 

The widths of the invariant mass spectra for simulation 
of W and Z production shown in Figure 4 are consistent 
with the inferred mass resolution from a study of inclusive 
two jet production at CDF ^^l- Over the dijet mass 
range from 60 to 126 GeV/c^ the fractional dijet mass 
resolution {5,n.- /trijj) is essentially constant at 15%. 



FITTING STANDARD MODEL PREDICTIONS 
TO THE DATA 

A sum of predicted background and resonant dijet 
mass distributions is fit to a histogram of the data in 11 
bins over the mass range 60 to 126 GeV/c^. The fitting 
procedure maximizes the likelihood function: 

L = exp(~xV2) 



where the x is given by 



1=1 



{N,~N,{f,b)f 



The optimal set of the two free parameters (/, 6) is 
found by minimizing the of the fit to the dijet mass 
distribution of the data. Ni and Ci are the number of 
events in the i*^ mass bin of the data, and the statistical 
error. The total number of events in the mass range from 
60-126 GeV/c^ is No = 2656. The data distribution 
dN/dmjj is described by a combination of background 
(B) and resonant W/Z signal (S). 



dN 

dm-ii 



= No{{l-f)B{m,,,b) + fS{m,,)) 
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Here, / is the fraction of the signal in data and b 
is the exponential slope of the background dijet mass 
distribution. Functions B{mjj,b) and S{mjj) represent 
parametrized shapes of the dijet mass spectrum for the 
background and the signal. The background distribution 
is represented by an exponential 

B{mjj~b) = Boe^"« 

while the signal is represented by the double Gaussian 

described above. The func- 
are normalized so that the 



distribution function S{mjj 
tions B{mjj,b) and S{mjj) 
area under each, over the mass range from 60 to 126 
GeV/c^, is one. 



126 



B{mjj, b)dm 



60 



■30 



S{mjj)dmjj = 1 



60 



The prediction for the number of events in the ith mass 
bin is: 



dN , 

dmjn 

dm,, ■ 



bin i "'"-JJ 



The cross section for production of and Zj {ow/z-i) 
can be expressed in terms of the fraction / 



Nof 



e 



WjZ-i' 



where L is the total integrated luminosity of 90 ± 4 pb^^ 
and Sv/jz-^ is the combined WjZ^ event detection 
acceptance — 0.023 ± 0.004. This final acceptance is an 
average of £w~i and £z^ quoted in Tables 1 and 2. The 
errors on the final acceptance are calculated as weigthed 
average of the errors on the individual 14^7 and Z7 ac- 
ceptances. 

From the LO QCD predictions for the background and 
signal cross sections, we would expect about 1% of the 
events in the fit region of the dijet invariant mass spec- 
trum (60-126 GeV/c^) to originate from the hadronic 
decay of W jZ bosons. The fit that maximizes the likeli- 
hood function gives / of — 0.05±0.05, corresponding 
to an unphysical negative cross section. This fit to the 
data is shown as the solid line in Figure 3. 

Since no signal is observed, we calculate upper limits 
on the W jZ^ cross section, cr, using a Bayesian proce- 
dure. After maximizing it with respect to &, the like- 
lihood function depends on the parameters cr, £, and 
evK/Z7- A posterior probability density is obtained by 
multiplying this reduced likelihood with truncated Gaus- 
sian prior densities for £ and twjz-^^ with mean and 
width equal to the central value and uncertainty of these 



parameters, respectively. The cross section is assigned a 
uniform prior. The posterior density is then integrated 
over the parameters L and jz-^^ a-nd the upper limit for 
the WjZ^ production cross section is obtained by calcu- 
lating the 95*'' percentile of this distribution (the value 
of cr such that 95% of the area under this distribution 
is below this value). The posterior density integrated 
over £ and Cw/z-^ as a function of cr is shown in Fig- 
ure 6. The result is a{pp — > W^) x BR{W jj) + 
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FIG. 6; Posterior density (Likelihood function multiplied by 
Gaussian priors for C and ew/z-y) integrated over £ and 
^w/z-y, versus a. 

a{pp Z-f) X BR{Z jj) < 50 pb compared to a 
Standard Model expectation of 13 pb. 



SUMMARY 

We have searched for the process pp — > W/Z + j + X 
in which the W and Z bosons decay hadronically. The 
sensitivity of the search is enhanced by fitting the dijet 
mass spectrum from pp 2 jet +7 to that expected 
from QCD models for the background and the line shape 
from the W/Z —> 2 jets decay The Standard Model 
prediction of the cross section for pp —>■ W/Z + 7, with 
the W/Z bosons decaying to quark- antiquark pairs, is 13 
pb. This cross section is for photons with of Et{"i) > 15 
GeV, Ai?(7 ~q)> 0.4 and for \/I > My - iTy. 

The data from 90 pb""'^ oipp interactions show no evi- 
dence of this process. By integrating the likelihood func- 
tion, we get a 95% confidence level upper limit to the 
cross section of 50 pb. This analysis method could prove 
effective for measuring the W/Z dijet mass signal in 
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this channel from a larger statistics sample of data being 
collected by the CDF experiment in Tevatron Run 2. 
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